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1.  INTRODUCTION 


The  application  of  Computational  Fluid  Dynamics  (CFD)  to  the  field  of  Projectile 
Aerodynamics  has  seen  continued  progress.  Early  Navier-Stokes  computational  techniques 
have  been  used  to  predict  pressure  distributions  on  projectile  configurations  at  transonic 
speeds  (Nietubicz,  Pulliam,  and  Steger  1980).  Later  studies  (Nietubicz  1981;  Nietubicz,  Inger, 
and  Danberg  1982;  Sahu,  Nietubicz,  and  Steger  1985)  have  provided  predictions  of 
aerodynamic  coefficient  data  for  conventional  projectiles,  hollow  projectiles  and  projectiles  with 
base  bleed.  Underlying  the  success  of  these  Navier-Stokes  solutions  was  the  initial 
generation  of  a  computational  grid. 

As  part  of  any  CFD  numerical  process,  the  region  in  which  the  equations  of  fluid  motion 
are  to  be  solved  must  be  discretized  into  a  finite  set  of  node  points.  At  each  of  these  nodes, 
the  laws  of  fluid  motion  are  enforced  with  an  appropriate  solver,  and  a  fluid  dynamic  solution 
can  be  obtained.  The  process  of  discretizing  this  flow-field  region  is  known  as  grid 
generation. 

Primary  methods  of  grid  generation  are  either  algebraic  or  are  based  on  solutions  of  partial 
differential  systems  of  elliptic,  parabolic,  or  hyperbolic  equations.  Work  on  projectile  grid 
generation  (Steger,  Nietubicz,  and  Heavey  1981)  included  both  algebraic  and  elliptic  grid 
generation  methods.  The  original  code  allowed  for  definition  of  projectile  shapes,  specification 
of  body  point  clustering  and  choices  of  grid  topology.  This  code  was  later  enhanced  to 
include  the  capability  to  generate  hyperbolic  grids  (Nietubicz,  Heavey,  and  Steger  1982).  The 
primary  advantage  of  the  hyperbolic  method  is  that  it  enforces  conditions  of  orthogonality  at 
each  grid  node  and  is  a  marching  procedure  as  opposed  to  a  more  time-consuming  elliptic 
method.  Although  this  increases  the  speed  and  flexibility  of  the  grid  generator,  the  impact  or 
need  for  this  method  was  not  fully  appreciated  until  computational  grids  were  required  for  the 
dome  bases  of  the  M864  and  M825  projectiles  (Nietubicz  and  Sahu  1988;  Sahu,  Nietubicz, 
and  Heavey  1988).  Up  to  this  point,  the  grid  generation  process  had  always  been  performed 
in  a  non-interactive,  batch  mode  environment. 

The  process  of  grid  generation  is  not  a  trivial  task  even  for  relatively  simple  shapes.  As 
body  configurations  become  more  complex,  the  task  of  creating  a  computational  grid  becomes 
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cumbersome  and  very  time  consuming.  The  grid  generation  process  previously  used  for 
projectile  aerodynamics  at  BRL  made  use  of  a  combination  of  batch  codes  and  static  graphics 
displays  which  provided  only  limited  analysis  capabilities. 

Through  the  introduction  of  powerful  graphics  work  stations  into  the  research  and 
development  workplace,  visual  analysis  has  become  an  important  analytical  tool  for  the 
engineer.  This  task  has  become  easier  to  perform  as  the  work  stations  have  become  faster  in 
their  compute  and  display  speeds.  Many  of  these  work  stations  support  zoom  and  pan 
capabilities  which  are  also  invaluable  in  the  analysis  process. 

Presently  there  are  several  general  purpose  grid  generation  codes  available  within  the 
United  States  which  can  handle  a  large  number  of  very  complex  configurations.  Many  of 
these  codes  are  well  known,  such  as  3DGRAPE  (Sorenson  1988,  1989),  EAGLE  (Thompson 
1988;  Thompson  and  Lijewski  1988),  GRIDGEN  (Steinbrenner,  Chawner,  and  Pouts  1989), 
GENIE  (Soni  1988),  and  VISUALGRID  (Cordova  1990).  These  codes  either  exist  already  in 
an  interactive  mode,  or  are  in  the  process  of  becoming  interactive. 

This  report  will  describe  the  evolution  of  the  BRL-PROGRID  (Projectile  Grid  Generation) 
code  to  an  interactive  platform,  and  will  discuss  some  of  the  new  enhancements  which  have 
been  added  to  the  code  for  projectile  grid  generation.  Additionally,  the  use  of  CAD  software  in 
providing  initial  body  point  distributions  will  be  discussed.  To  date  only  the  hyperbolic 
generator  has  been  fully  implemented  and  tested  on  both  projectile  and  store  separation 
configurations. 

2.  ORIGINAL  BRL-PROGRID  CODE 

2.1  History.  A  projectile  grid  generation  code  originally  developed  by  Joseph  Steger 
(Steger  1979)  was  written  in  FORTRAN  and  ran  on  a  VAX  VMS  system  in  a  non-interactive 
mode.  This  code  was  later  modified  to  include  a  hyperbolic  solver,  relaxation  parameter,  cell 
averaging,  and  secant  nose  configurations  (Steger,  Nietubicz,  and  Heavey  1981).  An  input 
data  file  was  read  and  a  data  file  containing  the  final  grid  was  produced.  Within  the  input  file 
the  user  could  select  from  an  algebraic,  elliptic,  or  hyperbolic  solver  to  generate  the  grid.  Also 
set  within  the  input  file  were  parameters  such  as  the  minimum  normal  spacing  at  the  body. 
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CDS;  spacing  at  the  outer  grid  boundary,  CDS2;  the  overall  normal  length  of  the  grid.  STOT; 
and  the  number  of  normal  points  in  the  grid,  KMAX.  If  any  of  these  parameters  required 
modification,  the  user  had  to  change  the  input  file  and  re-run  the  code.  Several  other 
parameters  are  also  set  in  the  input  file  depending  on  whether  the  body  is  being  input  as 
discrete  points  or  in  a  pseudo-analytical  representation. 

Next,  a  second  program  using  DISSPLA  software  graphics  calls  was  used  to  view  the  grid 
in  order  to  determine  if  it  was  acceptable  to  be  used  by  the  appropriate  solver.  This  process 
was  very  time  consuming  since  the  user  could  only  view  one  part  of  the  grid  at  a  time.  To 
complicate  this  task,  the  user  also  had  to  know  the  physical  coordinates  of  the  region  which 
was  to  be  viewed.  If  it  was  determined  that  the  grid  quality  was  not  acceptable,  the  input  data 
file  had  to  be  changed,  and  the  process  was  repeated. 

The  additional  requirement  for  grid  point  clustering  in  regions  of  interest  presented  another 
problem.  If  the  user  determined  that  it  was  necessary  to  cluster  grid  points  around  a  region  of 
rapidly  changing  geometry,  these  changes  had  to  be  made  in  the  input  data  file  and  it  became 
a  trial  and  error  iteration  until  the  desired  clustering  was  produced.  This  process  was 
repeated  for  each  zone  of  the  grid  until  the  entire  configuration  was  complete. 

Finally,  a  set  of  programs  was  used  to  create  the  grid  in  the  region  between  the  body  and 
the  base  and  then  combine  the  grid  zones  into  the  correct  format  for  different  solvers. 

Different  programs  were  used  depending  on  grid  zone  location,  desired  final  grid  orientation, 
and  the  solver  into  which  the  grid  would  ultimately  feed. 

2.2  Hyperbolic  Grid  Generation  Method.  For  most  projectile  applications,  the  outer 
boundary  is  unconstrained  and  simply  needs  to  be  placed  far  enough  away  from  the  projectile 
body  so  as  not  to  adversely  affect  the  flow  field  solution.  This  situation  represents  an  ideal 
application  for  a  hyperbolic  grid  generation  scheme. 

Once  the  body  points  have  been  adequately  defined  and  the  sting  or  cut  has  been 
determined,  a  grid  can  be  generated  using  a  hyperbolic  solver  similar  to  that  described  by 
Steger  and  Chaussee  (1980).  Before  the  actual  solver  can  be  implemented,  however,  the 
distance  to  the  outer  boundary  must  be  specified  and  either  constant  spacing  in  ti  or  some 
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type  of  stretching  function  is  required.  The  r\  stretching  used  here  is  determined  by  the 
following  relationship: 


As,  =  As,(1  =  . 


(1) 


Here  ASg  is  the  minimum  specified  grid  spacing  desired  at  the  wall  or  inner  boundary.  The 
parameter  e  is  determined  by  a  Newton-Raphson  iteration  process  so  that  the  sum  of  the 
above  increments  matches  the  known  arc  length  between  q  =  0  and  q  =  for  points  which 
have  the  same  value  of 

The  governing  equations  for  the  hyperbolic  solver  are  obtained  by  requiring  (1)  the 
coordinate  lines  ^  and  q  to  be  orthogonal  and  (2)  the  specification  of  a  cell  volume  or  area  for 
the  two  dimensional  case.  The  condition  of  orthogonality  requires 


A^  •  Aq  =  0.  (2) 

The  second  equation  is  obtained  by  specifying  a  grid  cell  volume  (or  area  in  two 
dimensions).  Since  the  grid  cell  volume  is  finite,  the  transformation  Jacobian  will  be  greater 
than  1,  i.e., 

dxdy  -  I  -  x^y^  \  d^d^  .  (3) 

The  set  of  grid  generation  equations  are,  therefore,  given  in  the  physical  plane  by 


=  0 


or  in  the  transformed  plane  by 


"  y^yn  =  0 

=  1/  V. 


(4) 


Using  local  linearization  for  this  set  of  nonlinear  differential  equations,  the  resulting  system 
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is  shown  to  be  hyperbolic  and  can,  therefore,  be  marched  in  the  direction. 


The  linearized  set  of  differential  equations  to  be  solved  numerically  is  written  in  vector  form 
as 


(5) 


where 


0  0 

.,0  ..0 

A  = 

.  B  = 

0  0 

0  0 
[yr, 

0 

1/  +  \/° 


and  x^,  y°,  etc.,  refer  to  known  conditions. 

The  set  of  Equations  (5)  are  solved  with  an  implicit  finite  difference  scheme  which  is  first 
order  accurate  in  the  r\  direction  (k)  and  where  central  differencing  is  used  in  the  ^  direction 
(y).  The  resulting  set  of  finite  difference  equations  becomes 


**i)  ^(0. 


*♦1. 


-  0>] 


2A^ 


Aq 


=  f. 


Rearranging  Equation  6  and  setting  Aq  =  A^  =  1.  results  in 


(6) 
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(7) 


^7^  =  r 

-2 


where 


{x^xo  .  yi“y°),. 


Equation  7  is  now  in  a  form  which  can  be  easily  solved  by  inverting  a  block  tridiagonal 
matrix  with  2x2  blocks.  The  terms  x°  and  y°  are  central  differenced  as 


^y.i, »  ■  ^j-i.k 
2 


//♦I.  A  “  Xy-i.  * 
2 


(8) 


The  terms  and  are  obtained  from  Equation  4  evaluated  at  the  old  station  (o).  That  is, 


0  0  0  0  r\ 

^  y^r  “8 

0  0  0  0  I  /  n 

Xn  -x„y^  =  V°. 


Solving  for  and  y^  with  and  y^  given  in  Equation  8  yields 


(9) 


(10) 
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Now  the  cell  volume  remains  to  be  specified.  This  specification  is  important  since  it  has 
the  effect  of  controlling  the  grid  evolution  as  the  solution  is  being  marched  out  from  the  body. 
The  method  chosen  here  is  straightforward  and  uses  the  stretching  function  given  by 
Equation  1 .  Specifying  the  minimum  spacing  at  the  wall  As^  and  the  total  number  of  points 
in  the  q  direction,  an  array  of  arc  lengths  As*  is  determined.  Since  the  Ax  is  known  along  the  j 
line,  the  volumes  are  calculated  by 


V  =  (x,.,.*-x,*).  (11) 

This  specification  of  cell  volumes  yields  smoothly  varying  grids  in  the  q  direction.  Grid 
volume  control  is  obtained  by  varying  the  arc  length  distribution  As*  and/or  surface  point 
distribution.  A  grid  generated  using  this  technique  is  shown  in  Figure  1.  Figures  2  and  3 
expand  the  base  region  and  show  the  grid  generated  for  the  concave  base  cavity  of  the  M864 
projectile.  An  additional  volume  specification  approach  can  be  found  in  Steger  and  Chaussee 
(1980). 

Due  to  the  concavity  of  the  base  region  on  the  M864  projectile  shown  in  Figure  3,  the  use 
of  a  hyperbolic  solver  would  normally  lead  to  a  grid  which  was  unacceptable  for  use  in  a  CFD 
code  solver.  The  grid  would  be  unacceptable  because  enforcing  conditions  of  orthogonality  in 
this  region  would  lead  to  grid  cells  crossing  and  overlapping  each  other.  A  more  general  form 
of  the  hyperbolic  equation  set  was  developed  by  Kinsey  and  Barth  (July  1 984)  which  enables 
grids  to  be  generated  for  complex  configurations  such  as  the  base  region  of  the  M864 
projectile.  The  generalization  is  obtained  by  replacing  the  backward  integration  in  q  of 
equation  6  with 


o/  \ 

~  0.*) 

Aq 


(1  -aK'',)* 


(12) 


where  a  is  used  to  control  the  type  of  finite  difference  marching  algorithm  and  has  the 
resultant  effect  of  relaxing  the  degree  of  orthogonality  of  the  grid.  When  a  =  1 ,  the  original 
scheme  is  maintained  and  the  grid  is  everywhere  orthogonal.  For  a  >  1  the  numerical  error 
term  is  dissipative  and  the  grid  is  no  longer  fully  orthogonal. 
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The  final  grid  generated  using  the  generalized  equations  is  shown  in  Figure  3  and  required 
ALPHA  set  at  6.  The  effect  of  lowering  the  value  of  ALPHA  to  3  is  shown  in  Figure  4.  Any 
further  attempt  to  reduce  the  value  of  ALPHA  leads  to  grid  line  crossing  for  this  configuration. 

An  additional  change  was  made  in  the  code  to  reduce  the  propagation  of  body 
discontinuities  into  the  grid  as  the  grid  is  marched  toward  the  outer  boundary.  This  is 
accomplished  through  a  smoothing  parameter,  ITERV,  which  is  used  to  average  the  grid  cell 
volumes  of  Equation  1 1 .  ITERV  is  the  number  of  averaging  sweeps  taken  for  the  specified 
volume.  Wj .  An  ITERV  value  of  36  was  required  to  generate  a  smooth  grid  with  only  minor 
propagation  of  the  concave  region  of  the  base  surface  toward  the  outer  boundary  of  the  grid 
(see  Figure  2).  Setting  ITERV  equal  to  0  disables  this  averaging  procedure  and  generates  the 
grid  shown  in  Figure  5.  Here,  the  effect  of  the  base  cavity  is  seen  to  propagate  toward  the 
down  stream  boundary. 

3.  INTERACTIVE  BRL-PROGRID  CODE 

3.1  Body  Definition.  One  of  the  first  difficulties  encountered  in  trying  to  develop  a  finite 
difference  grid  is  that  of  attempting  to  get  an  accurate  definition  of  the  body  surface.  In  an 
effort  to  solve  this  problem,  it  was  decided  that  a  CAD/CAM  program  could  be  used  in 
conjunction  with  the  engineering  drawings  to  produce  an  accurate  representation  of  the 
configuration  of  interest.  A  PC  program  called  CADKEY  was  chosen  since  it  was  easy  to  use 
and  provided  output  files  in  a  format  which  could  easily  be  used  to  create  a  discrete  set  of 
body  points. 

Engineering  drawings  are  rapidly  and  easily  entered  into  the  CADKEY  program  via  a 
combination  of  mouse  and/or  keyboard  input.  An  example  of  a  drawing  created  using 
CADKEY  is  shown  in  Figure  6.  Once  the  design  is  inside  CADKEY  it  can  be  written  out  in  a 
variety  of  different  formats.  One  of  these  formats  is  a  ".CDL"  file.  This  file  contains  the 
analog  information  defining  each  of  the  curves  which  make  up  the  projectile’s  surface.  A 
sample  ".CDL"  file  for  the  M864  projectile  configuration  is  shown  in  Figure  7.  A  FORTRAN 
program  is  then  run  which  takes  the  information  contained  in  the  ".CDL"  output  file  and 
calculates  a  set  of  discrete  points  to  define  the  projectile  surface.  The  data  file  containing 
these  points  is  then  used  as  part  of  the  input  file  to  the  interactive  grid  generator. 
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3.2  Grid  Generation.  The  process  of  continually  entering  and  exiting  several  different 
programs  in  order  to  make  even  relatively  simple  changes  to  the  generated  grid  was  very  time 
consuming  and  inefficient.  In  order  to  speed  up  this  necessary  and  important  process, 
interactive  graphics  were  incorporated  into  the  grid  solver.  The  developed  graphical  interface 
makes  use  of  the  Silicon  Graphics  GL  Library  and  the  geometry  engines  of  the  Iris  4D  series 
work  stations  to  provide  more  rapid  computation  and  analysis  of  a  hyperbolic  grid. 

Incorporated  into  the  interactive  grid  generator  was  the  ability  to  zoom  and  pan  the  final 
grid.  This  allows  the  user  to  rapidly  view  a  region  of  the  grid  which  might  be  suspected  of 
having  discontinuities.  The  user  can  also  separately  control  the  sensitivity  of  the  mouse  for 
zooming  and  panning.  Color  was  also  added  to  facilitate  the  analysis  of  zonal  grids.  The  grid 
shown  in  Figure  1  consists  of  three  zones  when  it  is  originally  generated:  the  body  grid, 
shown  in  red;  the  base  grid,  shown  in  yellow;  and  the  "fill"  grid,  shown  in  blue.  Without  the 
use  of  color,  the  process  of  checking  for  correct  grid  overlap  and/or  meshing  would  be  a  very 
difficult  task.  Another  feature  which  was  added  was  the  ability  to  interactively  redistribute 
points  along  the  body  surface.  This  gives  the  user  unlimited  flexibility  to  add,  subtract,  or 
recluster  the  points  defining  the  projectile  surface  and  immediately  regenerate  a  new  grid. 

Also  added  was  the  ability  to  interactively  change  the  parameters  CDS,  CDS2,  STOT,  ALPHA, 
ITERV,  and  KMAX  which  provide  conditions  or  constraints  for  the  formulation  of  the  grid. 

Once  any  or  all  of  these  parameters  have  been  changed,  the  user  can  immediately 
regenerate  the  grid  and  view  it  to  see  if  the  changes  brought  about  the  desired  results.  The 
user  can  also  create  the  grid  region  between  the  body  and  the  base  with  one  keystroke  once 
the  body  and  base  grids  have  been  created  and  saved.  These  new  abilities  provide  a  more 
efficient  method  of  creating  and  analyzing  a  computational  grid.  Every  part  of  the  generation 
process  is  handled  by  one  program,  rather  than  using  several  programs  as  was  previously 
required.  This  makes  it  easier  for  the  user  to  more  rapidly  create  a  finite  difference  grid. 

The  opening  screen  for  BRL-PROGRID  is  shown  in  Figure  8.  The  Main  Menu  options 
currently  available  are:  Input  Body  File,  Input  Base  File,  Re-Distribute  Points,  Change 
Parameters,  Form  Grid,  Un-Do  a  Regrid,  Show  Body  Only,  Show  Base  Only,  Show  Body  and 
Base,  Save  Grid,  Exit  Menu,  and  Exit  Program.  A  more  detailed  discussion  of  these  menu 
options  and  a  sample  session  are  presented  in  Appendix  A  which  also  serves  as  a  Users’ 
Guide  for  the  program. 
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4.  GRID  POST-PROCESSING 


Once  the  user  has  created  and  saved  all  the  desired  grid  zones,  BRL-PROGRID  is  closed 
and  BRL-GRIDCOMB  is  run.  BRL-GRIDCOMB  is  written  in  FORTRAN  on  an  Iris  4D  work 
station,  and  uses  the  Iris  GL  Library  and  geometry  engines  for  controlling  the  graphics.  This 
program  takes  as  input,  multiple  PLOT30  (Walatka  and  Buning  1989)  format,  2D,  single  grids 
and  allows  the  user  to  interactively  combine  all  the  grids  into  the  correct  format  and  orientation 
for  either  the  F3D  code  (Sahu  1988),  or  the  AXISYMMETRIC  Navier-Stokes  codes  (Nietubicz, 
Pulliam,  and  Steger  1980).  It  also  allows  the  user  to  interactively  split  up  a  large  grid  into 
smaller  grids  for  input  into  the  zonal  F3D  code.  When  the  user  is  done,  the  code  writes  out 
the  final  grid(s)  in  the  desired  format.  BRL-GRIDCOMB  also  allows  the  user  to  zoom  and  pan 
across  the  grid(s). 

5.  SUMMARY 

An  interactive  hyperbolic  projectile  grid  generation  capability  has  been  developed  for  use 
on  the  Iris  Silicon  Graphics  workstations.  The  addition  of  user  interaction  through  graphics 
has  made  the  hyperbolic  grid  generation  task  easier  and  faster  to  complete.  An  existing  batch 
code  was  ported  to  the  Iris  Silicon  Graphics  4D/xxx  series  of  work  stations  in  order  to  facilitate 
the  incorporation  of  interactive  graphical  analysis  capabilities  into  the  grid  generation  process. 
Capabilities  gained  from  this  port  include  zoom/pan  ability,  body  surface  point  redistribution, 
and  spacing  control  in  the  normal  direction.  All  of  these  tasks  can  now  be  completed  without 
the  user  being  required  to  exit  the  program,  make  changes  in  the  input  file,  and  enter  the 
program  again  with  the  changed  input  file. 
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APPENDIX: 

BRL-PROGRID  USER’S  GUIDE  FOR 
INTERACTIVE  HYPERBOLIC  GRID  GENERATION  PROGRAM 
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Intentionally  left  blank. 
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INTRODUCTION 


1)  Type  "progrid"  at  the  command  prompt. 

2)  The  Main  Menu  will  appear  in  the  upper  right  corner,  click  the  right  mouse  button  on  "input 
body  file." 

3)  Control  will  move  to  the  textport  at  the  bottom  of  the  screen  where  you  can  type  in  the 
name  of  your  input  file.  (The  structure  of  the  input  file  is  outlined  later  in  the  Input  File  Format 
section.) 

4)  If  the  data  file  is  read  in  correctly,  the  body  and  grid  will  appear  in  the  viewport. 

5)  Use  the  right  mouse  button  to  translate  the  grid,  the  middle  mouse  button  to  zoom  the  grid, 
and  the  left  mouse  button  to  reset  the  translation  and  zooming  to  default. 

6)  If  the  computed  grid  is  unsatisfactory,  press  the  "m"  key  to  activate  the  Main  Menu  and 
perform  one  of  the  following  options; 

a)  Click  the  right  mouse  button  on  the  "Change  Parameters"  option.  Control  will 
transfer  to  the  textport.  The  parameters  you  can  change  are: 

cds  -  spacing  at  body  (-1  for  constant  distance) 

cds2  -  spacing  at  outer  body 

stot  -  total  distance  to  outer  boundary 

alpha  -  smoothing  parameter  to  "relax"  grid  orthogonality  conditions 

iterv  -  number  of  grid  smoothing  iterations 

kmax  -  number  of  grid  points  in  the  normal  direction. 

If  you  do  not  wish  to  change  the  listed  parameter,  hit  the  "return"  key,  and  the  old 
value  will  remain  unchanged.  To  enter  a  value  of  zero  type  in  (-9999.0).  Changes  will  be 
saved  in  a  file  called  "new_param." 

b)  Click  the  right  mouse  button  on  the  "Re-Distribute  Points"  option.  Refer  to  the 
sample  point  redistribution  session  later  in  this  text  for  more  information. 

7)  Return  control  to  the  Main  Menu  and  click  on  the  "Form  Grid"  option  to  compute  a  new 
grid  using  the  new  parameters  or  new  point  distribution. 

8)  If  the  grid  is  satisfactory,  turn  on  the  Main  Menu  (MKey)  and  click  on  the  "Save  Grid" 
option.  Control  will  pass  to  the  textport  where  you  should  enter  a  1  to  save  your  new  body 
grid. 

9)  The  grid  is  created  in  two  dimensions,  however,  you  can  save  the  grid  in  either  ASCII  or 
"unformatted"  PLOT3D  -  2D  or  3D  format.  Also,  at  this  point,  the  grid  can  be  rotated  into  a 
3D  solid;  first  enter  the  total  number  of  degrees  of  rotation,  then  the  number  of  planes 
desired.  The  rotated  output  file  can  be  saved  in  X-Y-Z  columnar  format,  and/or  PLOT3D 
format. 
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1 0)  Click  on  “Exit  Program"  and  confirm  in  the  textport  to  exit. 


BRL-PROGRID  INPUT  BODY  FILES  and  DATA  FILE  STRUCTURE 

1)  The  input  file  structure  of  the  BRL-PROGRID  program  has  been  changed  so  as  to  keep 
compatibility  with  the  VAX  batch  code  BRL  grid  generator. 

BRL-PROGRID  input  file  structure: 


variable  or  flag  name  data  type 
isolv,  i3d,  ibase,  nd,  Imax,  alpha,  iterv 

jmax,  kmax,  ibc,  idat,  iscal,  icius.  stot,  cds,  cds2 

xx(i)  yy(i) 

Here,  isolv,  i3d,  ibase,  nd,  Imax,  alpha,  iterv,  jmax,  kmax,  ibc,  idat,  iscal,  and  icius  are 
integers;  and  stot,  cds,  cds2,  xx(i),  and  yy(i)  are  reals. 

These  inputs  have  the  following  meanings: 

isolv  -  0  for  elliptic  solver  (currently  an  unsupported  feature) 

- 1  for  hyperbolic  solver  (always  use  integer  1) 

i3d  -  0  for  2D  grid  (no  spin)  -  most  commonly  used  setting 

-  1  for  3D  grid  (half  plane  spin) 

-  2  for  3D  grid  (full  plane  spin) 

ibase  -  0  for  no  base  grid  (most  used  option) 

-  1  for  exponential  clustering  in  base  region 

-  2  for  Vinokur’s  clustering  in  base  region 

nd  -  number  of  circumferential  planes  (usually  zero) 

Imax  -  number  of  circumferential  planes  (same  as  nd) 
alpha  -  smoothing  parameter  to  "relax"  orthogonality  conditions 
iterv  -  number  of  grid  smoothing  iterations 
jmax  -  number  of  longitudinal  planes 
kmax  -  number  of  normal  planes 

ibc  -  0  for  boundary  conditions  along  negative  x-axis  at  J=1 

-  1  for  B.C.'s  perpendicular  to  x-axis  at  J=1 

-  2  for  boundary  conditions  along  positive  x-axis  at  J=1 
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idat  -  0  for  body  points  created  analytically  within  program  (rarely  used  setting) 

-  1  for  body  points  read  in  from  input  file  (unit  5)  -  most  commonly  used  setting 

-  2  for  body  points  read  in  from  separate  input  file  (unit  7) 

iscal  -  0  for  no  rescaling  of  hyperbolic  grid 

-  1  for  rescaling  of  hyperbolic  grid 

icius  -  0  for  original  clustering  in  normal  direction 

-  1  for  Vinokur’s  clustering  in  normal  direction 

stot  -  total  distance  to  outer  boundary  (-1  for  constant  spacing  using  (kmax  *  cds)  to 
determine  the  final  distance  to  the  outer  boundary) 

cds  -  spacing  at  body 

cds2  -  spacing  at  outer  boundary 

xx(i)  -  X  position  of  body  surface  points  (j  =  1  to  jmax) 

yy(i)  -  Y  position  of  body  surface  points  0  =  1  to  jmax) 

(Note:  all  read  statements  are  free  formatted.) 


BRL-PROGRID  MAIN  MENU  PANEL 

Input  Body  File:  Used  to  input  a  body  shape  from  an  existing  file.  The  user  can  type  "quit" 
here  to  quit  the  program. 

Input  Base  File:  Used  to  input  a  base  shape  from  an  existing  file.  (A  base  is  a  body  rotated 
90  degrees.)  The  user  can  type  "quit"  here  to  quit  the  program. 

Re-Distribute  Points:  Used  to  change  the  point  distribution  of  the  body  or  base  file.  (An 
example  of  how  to  use  this  option  follows  this  section.) 

Change  Parameters:  Changes  the  input  parameters:  cds,  cds2,  stot,  alpha,  iterv,  and  kmax. 

Form  Grid:  Forms  new  grid  after  changing  an  input  parameter,  or  redistributing  the  body  or 
base  points. 

Un-Do  a  Regrid:  Un-Does  a  regrid. 

Show  Body  Only:  Displays  only  the  BODY  (red  grid  in  Figure  1)  for  manipulation  if  both  a 
body  and  base  grid  have  been  created. 

Show  Base  Only:  Displays  only  the  BASE  (yellow  grid  in  Figure  1)  for  manipulation  if  both  a 
body  and  base  grid  have  been  created. 
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Show  Body  and  Base:  Show  both  the  BODY  and  BASE  if  a  body  and  base  grid  have  been 
created. 

Note:  If  you  select  this  option,  it  is  only  possible  to  zoom  and  translate  the  grid. 
Parameter  changes  and  point  re-distribution  are  not  possible. 

Save  Grid:  Saves  the  grid  as  a  PLOT3D  format.  2D  or  3D,  ASCII  or  unformatted  data  file. 
Also  contains  an  option  for  3D  grid  rotation. 

Exit  Menu:  Turns  the  Main  Menu  off,  and  sends  mouse  control  to  the  viewscreen.  (The  Main 
Menu  must  be  turned  off  to  zoom  and  translate.) 

Exit  Program:  Clears  the  textport  and  exits  the  program  if  you  confirm  your  selection  in  the 
textport. 


MOUSE  BUTTON  FUNCTIONS 


Left  Mouse  Button: 
Right  Mouse  Button: 
Middle  Mouse  Button: 


Reset  zoom  and  translate  to  default. 
Translate  object. 

Zoom  object. 


BRL-PROGRID  KEYBOARD  MACROS 

(These  keys  are  activated  when  the  Main  Menu  is  turned  off.) 

Escape  key:  Exits  program. 

B  key:  Turns  the  viewport  background  color  black  (default). 

C  key:  Displays  the  current  X  and  Y  translations  and  current  zoom  factor  in  the  textport. 

F  Key:  Creates  a  "fill  grid"  (blue  grid  in  Figure  1)  between  a  body  grid  and  base  grid. 

H  key:  (Help  Key)  Draws  a  picture  of  the  mouse,  telling  how  to  zoom,  translate  and  reset 
transformations. 

L  key:  Gives  the  current  X  and  Y  translations  and  zoom  factors,  and  asks  for  new  values. 
This  option  will  allow  keyboard  or  data  file  entry  of  translation  and  zooming  factors 
without  using  the  mouse.  If  you  do  not  wish  to  change  a  value,  hit  return,  and  the  old 
value  will  remain,  if  you  do  change  a  value,  make  sure  you  enter  it  as  a  real  number 
(include  decimal  point).  For  data  file  input,  the  user  is  prompted  for  an  input  file 
name.  The  format  for  this  file  is:  x-translation,  y-translation,  zoom  factor  in  free- 
format. 

M  key:  Turns  the  Main  Menu  panel  on. 

Note:  Menu  panel  must  be  turned  off  in  order  to  zoom  and  translate  the  grid. 
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P  Key:  Position  set/toggle.  Using  this  key,  the  user  can  set  two  viewing  locations  and  toggle 
between  them.  The  user  is  prompted  for  input  as  to  whether  they  are  setting  a  toggle 
position  or  are  toggling  between  two  positions  which  are  already  set. 

R  key:  Toggles  grid  reflection  on  and  off. 

S  key:  Save  a  zoomAranslation  data  file. 

W  key:  Turns  the  viewport  background  color  white. 

(up  arrow):  Increase  zoom  sensitivity. 

(down  arrow):  Decrease  zoom  sensitivity. 

(right  arrow):  Increase  translate  sensitivity. 

(left  arrow):  Decrease  translate  sensitivity. 


Example  of  Point  Redistribution  Menu  Option: 

1)  Click  the  right  mouse  button  on  the  "Re-distribute  Points"  option  in  the  Main  Menu  after 
reading  in  a  body  or  base  file. 

2)  Click  the  left  or  right  mouse  buttons  to  move  along  the  body’s  points.  Move  to  the  first 
desired  point  and  click  the  middle  mouse  button.  The  selected  point  will  turn  yellow  to 
indicate  that  it  has  been  selected  as  an  endpoint.  Now  move  to  the  next  desired  point  and 
click  the  middle  mouse  button  to  select  this  point  as  your  second  endpoint. 

Note:  You  can  only  select  two  (2)  endpoints  to  re-distribute  between.  To  un-select  an 
endpoint,  move  to  that  point  and  click  the  middle  mouse  button.  The  point  will  no  longer  be 
displayed  in  yellow,  and  you  can  now  select  a  new  endpoint. 

3)  Hit  the  R-Key  to  redistribute  points  between  the  chosen  endpoints.  You  will  be  prompted 
for  input  as  to  how  many  points  you  would  like  to  have  between  the  chosen  endpoints. 

4)  You  will  now  be  asked  if  you  want  linear  or  Vinokur’s  clustering  in  the  X-direction.  If  you 
choose  Vinokur’s,  you  will  have  to  input  "dyO"  and  "dyl",  which  are  user  specified  endpoint 
spacings  at  the  left  and  right  endpoints  respectively. 

5)  If  this  new  point  distribution  is  satisfactory,  type  the  K-Key  to  keep  the  new  points,  and  hit 
the  Escape-Key  to  go  back  to  the  Main  Menu.  Click  on  "Form  Grid"  to  compute  the  new  grid. 
If  you  wish  to  un-do  the  redistribution,  hit  the  U-Key  and  either,  re-do  the  redistribution  or 
press  the  Escape-Key  to  go  back  to  the  Main  Menu. 
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Re-Distribute  Points;  key  functionality: 

Left  Mouse:  Decreases  the  point  number  index. 

Right  Mouse:  Increases  the  point  number  index. 

Middle  Mouse:  Selects/Unselects  the  current  point  as  an  endpoint. 

K-Key:  Keep  the  point  redistribution. 

R-Key:  Redistribute  the  points. 

S-Key:  Show  currently  selected  redistribution  end  points. 

U-Key:  Undo  the  point  redistribution. 

up-arrow:  Places  the  user  in  "fast  scroll"  mode  for  left  and  right  mouse  buttons.  In 

this  mode,  the  user  can  "press-and-hold"  a  mouse  button  for  rapid 
increase  or  decrease  of  point  number  index. 

down-arrow:  Places  the  user  in  "slow  scroll"  mode  for  left  and  right  mouse  buttons. 

In  this  mode,  the  user  must  press-and-release  the  mouse  button  each 
time  he  wishes  to  increase  or  decrease  the  point  number  index. 

Escape-Key:  Exits  the  Re-distribute  subroutine,  and  goes  back  to  the  Main  Menu. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it 
publishes.  Your  comments/answers  below  will  aid  us  in  our  efforts. 

1.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of 
interest  for  which  the  report  will  be  used.)  _ _ 


2.  How,  specifically,  is  the  report  being  used?  (Information  source,  design  data,  procedure, 
source  of  ideas,  etc.)  _ _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or 
dollars  saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please 
elaborate.  _ _ 


4.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports? 
(Indicate  changes  to  organization,  technical  content,  format,  etc.)  _ 
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